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ABSTRACT 


This work reports theoretical investigations of the surface motion 
in solids caused by reflection of an intense acoustic pulse at a pressure. 
release surface done using the surface velocity reflection coefficient. 
A formula was derived for the surface velocity reflection coefficient 
which relates the surface particle velocity for an incident compression 
wave to the particle velocities of the incident and reflected waves. 

To obtain the surface profiles, two different assumptions were used. 
The first assumption was the conventional one. In this assumption, 
Surface motion was neglected during the period that the pulse strikes 
the surface and it is called the fixed surface approximation. The 
second assumption is to consider the motion and distortions of the 
Surface and is called the moving surface approximation. 

Four different materials were investigated and the results were 


compared for different materials and in the two different approximations. 
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I. INTRODUCTION 


Impulsive loads such as are produced by explosions and impacts may 
cause the failure of materials. An impulsive load is distinguishable 
from a conventional load primarily by its suddenness of application and 
its brevity of duration. Because of these loading characteristics, the 
reaction of a solid body to an impulsive load must be understood in 
terms of the propagation of waves in the body generated by the loading. 
If tne loading is of sufficiently large amplitude, these waves produce 
extensive fracturing and large permanent distortions in the body upon 
which it acts. The phenomena associated with failure under impulsive 
loads are markedly different from those connected with failure under 
conventional loads. The reaction of materials to the extremely high, 
very transitory pressures imposed on them by impulsive loading is of 
great fundamental importance [Ref. 1]. 

The basic character of the fracture and plastic flow which are 
caused by impulsive loading is beginning to be understood [Ref. 1]. 

A comnon fracture process results by the generation of shear waves 

whose amplitude exceeds the shear strength of the material in which 

the wave is propagating. Sucn shear waves are usually generated by 

the reflection of a strong compressive wave from boundaries. Anotner 
common fracture process results when the pressure release wave generated 
by the reflection of the compressive wave from a boundary has too great 
an amplitude. A third common failure process involves the superposition 
of shocks reflected from various boundaries. Thus, understanding the 
way in which intense waves damage solids requires an understanding of 


reflection of waves from boundaries. 
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In the limited theoretical discussion which has been given such 
processes to date [Ref. 5], the reflection process is discussed using 
expression obtained for sivall amplitude waves. Experiments show that 
this discussion is sufficient for qualitative, but not always quantita- — 
tive, understanding of some, but not all, of the observed phenomena. 

The laws governing oblique reflection of elastic waves at a free 
surface are well known [Refs. 2, 3]. There is a particle motion 
associated with the reflection of these elastic waves that strike 
free surfaces obliquely. Tne object of this thesis is to study 
possible departures from the behavior predicted by the reflection 
formulae which result because of the finite amplitudes of waves 
involved by investigating the surface motion in solids caused by 
reflection of an intense acoustic pulse at a pressure release 


surface. 











IT. SURFACE ii0TION DUE TO REFLECTION OF ELASTIC WAVES 


A. ELASTICITY 

Mathematical theories of elasticity are developed in the same 
general way. First, the notations of stress and strain are developed. 
Second, a stress-strain relationship between these quantities, or 
their derivatives, is assumed which idealizes the behavior of actual 
materials. Finally, using this relationship, equations of motion or 
equilibrium are set up which enable the state of stress or strain to 
be calculated when a body is subject to prescribed forces. 

In general, the stress on a surface element in a solid body 
does not act normally to the surface. It has components both normal to 
the plane and tangential to it. In Fig. 2-1, the coordinate system has 
three mutually perpendicular axes Cvs Oy» 0, The stresses acting on 
three planes normal to these axes which pass through a point P are con- 
Sidered. There will be nine components of stress. These are denoted 
, etc. The first letter in the suffix indicates 


by o 


» O 


memayy? 922° Xz 


the direction of the stress and the second one defines the plane 

in which stress is acting, If an infinitesimal rectangular parallelepiped 
around P, with its faces normal to the axes, is considered, by taking 
moments, 1t may be seen that for equilibrium o, =o6,., 06._ = o_., and 


xy Ye XZ ZX 
Typ Toys Thus there are only six independent components of stress. 
These completely define the stress at the point. 

The strain of the medium is the relative displacement of the material 
particles in the medium. A rigid translation or rotation of a medium 
is not considered a strain. In order to find the strain at the point, 


its change in position relative to adjacent points must be analyzed. 
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Figure 2.1. Stress components acting on an infinitesinal rectangular 
parallelepiped 
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Consider a point very close to P, which before displacement had 
coordinates (x t dx), (y t dy), (z+dz). Let the displacement which 
it has undergone nave components (u + éu, v + 6v, wt dw), where u, Vv 
and w are the components of the displacement of point P. If &x, dy 


and 6z are sufficiently small. 


dU dU 


iain 
_ av av 
(2-1) WAS coe Oh! - weap s 62 


Sy = 2H gx 4 2 vate $z 
aX Jy 
The relative displacement of all surrounding points P may be found 


; oz... OU OU OU OV BV dV dw dw dw 
if the values of the nine quantities IX? By? BZ? BX? ay? OZ? OX? ay? OZ 


at the point are known. These nine quantities are regrouped as follows: 


ql oul! = eee = 2 
XX 9X yy ay red wee 
ow dV dU CW dV dU 
™ ee ee = —+ — =——t — 
(2-2) Arius EN Ope Teh Cyl Oye “xy Ox. OY 
— 2 OV oe 1 aul mw og LBV | BU 
eu, = oy” 2z Cu, = 37 7 3x cw, = ox 7 ay 
The first three quantities, Ex? Ey? €2 7 correspond to the fractional 
expansion. The second three my =e and Ey correspond to the components 


of shear strain in the planes denoted by their suffixes. The last three 
Wy» w and a, do not correspond to a deformation of the element around P, 
but are the components of its rotation as a rigid body. 

Experimentally it is found that for most solids, the observed strains 
are proportional to the applied load, provided that the load does not 


exceed the elastic limit. This law is known as generalized form of 
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Hooke's law, and can be stated as: Each of the six components of stress 
is at any point a linear function of the six components of strain. This 
law may be written as 


Cratyz * (15Ezx * —16Exy 
a 


Oxy = CrExx * Cy2eyy * ©13&22 


“ + 


Ovy = CorExyx 7 Cantyy * Ca3&27 * Coat yz * C25Ezx * ©26&xy 


| + 


Foz ~ S31Exx * S32&yy * ©33&22 * S3qeyz * ©35Ezx * °36°xy 
(2-3) 
Ov ~ SarExx * Cqztyy * 43®2z * Cqgeyz * S45 2x * “aE xy 


+ + 4 + 


Fox = S1Exx * S52kyy * ©5322 * S5afyz * S55f2x * —56Exy 


€ € € e. 
63 zz y “64 yz * ©65 2x * “66 xy 


= € € 
Syy ~ 61 xx ey, yy 
The coefficients in these relations are known as the elastic constants 
of the material. 


In an isotropic solid, the values of the coefficients are independent 
of the set of rectangular axes chosen. For such a solid, there are only two 


independent elastic constants, » and pw known as Lamé's constants. All the 
elastic constants can be expressed in terms of the Lamé's constants. 
Thus, the Lamé's constants completely define the elastic behavior of an 


isotropic solid. For such a material, equation 2-3 may be written 


0 


+ = + = + 
o AA Cue o AA 2ue Oo AA eve, 


yy yy 22 
(Oy, 


0 HE 0 HE 0 = HE 


yz yz ZX ZX xy xy 
Here A = ene 4 Evy i E5y and is called dilatation. It represents the 
change in volume of a unit cube. 

For convenience, four elastic constants are normally introduced. 


One is Young's modulus E, which is given by Cue leKy so that: 


i 





Mees) og = BASLE Zu) 


A +p 


A second is Poisson's ratio v defined as era eo so that: 


(2-6) v= 


A third is the bulk modulus k defined as the ratio between the applied 
pressure and the fractional change in volume when the solid is subjected 


to uniform hydrostatic compression so that: 


eu 


(2-7) SUS ae Ge 


Finally, the rigidity of shear modulus, which is identical with Lamé's 
constant u, is defined to be the ratio between the shear stress and the 


Shear strain. 


B. EQUATION OF MOTION IN AN ELASTIC MEDIUM 

In order to obtain the equations of motion. for an elastic medium, 
the variation in stress across a small parallelepiped with its sides 
parallel to a set of rectangular axes are considered (Fig. 2.2). The 
components of stress will vary across the faces; to Resa the force 
acting on each face, the value of stress at the center of each face 
times the area of the face has been taken. 

Six separate forces will act parallel to each axis. The resultant 
force acting in the x-direction is 


Wy 


+ — “ 
(cee a6 6x) 6y6Z - 0, dy 6z 





“7 (ore oe eo yy éy)6xéz = 6 5x82 + 


XY ay 
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Figure 2.2. Stresses acting on a smal] rectangular parallelepiped 
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By Newton's second law of motion, neglecting body forces such as 
megravity, this will be equal to 


2 
(pSx sysz) 24 


at? 


where p is the density of the element and u is the displacement in the 


x-direction. Thus the equations of motion are 











2 3 3 ,) 
(2-8) 52. = Oe Me ae 
2 OX oy OZ 
at 
Similarly 
2 oO. 00 30 
at aX ay Zz 
2 30 re) 30 _ 
(2-10) 9W._ 4X4 _ 24 4 ee 
st OX oy OZ 


C. ELASTIC ilAVES 
Equations 2-8, 2-9 and 2-10 will hold whatever the stress-strain 
behavior of the medium. By substituting the equations 2-4 for the 


Stress components in Equation 2-8, one obtains 


ae S (0A + 2yue on) + ay HEY + ae ie 


By using the equations 2-2, one finds 


2 
ou _ oA 2 
(2-11) ve a (A + Mee + uVou 
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Similarly 


2 
V 3A 2 

(2-12) Zz = (A+ lee + uv u 
=13) au (a+ poh + uv? 

= io ulay T WV w 


—, 
Equations 2-11 througn 2-13 are the equations of motion of an 


isotropic, elastic solid in which body forces are absent. They predict 
the propagation of two types of waves through the medium, If both 
Sides of the above equations are differentiated with respect to x, y and 


z respectively, they become 


2 
(2-14) oS = (A + 2u)v%a 
r) 


This is the wave equation and shows that tne dilatation A is propagated 


through the medium with velocity 
1 
(2-15) c, = L(A + 2u)/o]? 


On the other hand, if 4 is eliminated, it can be shown that 
ek 

(2-16) 5 = va, 

ot 


where wy is the rotation about the x-axis. Similar equations may be 


obtained for wy and Wa Thus the rotation is propagated with velocity 


(=17) Cy = (w/o) 


This shows that in the interior of any elastic solid waves may be 
propagated with two different velocities. These two types of waves 


should be termed dilatation waves and distortion waves. 
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Any plane wave propagated through on isotropic elastic medium 
must travel one or the other of the above velocities. The waves which 
propagate with velocity Cc, are called longitudinal waves. The motion 
is along tne direction of propagation. The other type waves, which 
propagate with velocity Co» are called transverse waves since the 
motion is parallel to the wave front. 


The equations 2-14 and 2-16 are of the form 


(2-18) 3¢ 


t 


cob) 


When the deformation is a function of only one coordinate, e.g. x, the 


equation beconss: 


2 
23 

(2-19) ary = ¢ 
ate ax 


The general solution of this is: 


= f(x - ct) + F(x + ct) 


where f and F are arbitrary functions depending on initial conditions, 
f corresponding to a plane wave travelling along the positive direction 


of the x-axis, and F to one in tne direction opposite to tnis. 


D. REFLECTION AT PRESSURE RELEASE SURFACE 

Normal incidence. Any elastic wave will be reflected when it 
reaches a free surface of the material in which it is travelling. The 
simplest case occurs when the wave strikes the surface normally. In 
a longitudinal wave, since the stress normal to the surface at the 


Surface must be zero, the reflected pulse must be opposite in sense 
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to the incident pulse (compression reflected as tension and vice 
versa). | 

Oblique incidence. When an elastic wave strikes a free surface 
obliquely, the situation is much more complicated. This complication 
arises from the fact that the energy of incident wave is partitioned 
into two reflected waves instead of only one. When either a longitudinal 
or a transverse pulse strikes a surface obliquely, both a longitudinal 
and a transverse wave are reflected. (Figure 2.3). 

Consider first a longitudinal wave striking a surface obliquely. 
There will be a reflected longitudinal portion in which the angle of 
incidence a is equal to the angle of reflection. There will also be 
a transverse wave which is generated at the point of reflection. 

From Snell's law, the relationship between the angle of incidence a 


and the angle of reflection 6 is 


P C 1 
(2-20) Sins > = [2(1 - v)/(1 = 2v) 7° 


From equation 2-20, it is seen that the angle of reflection g is the 
function of the incidence angle and tne Poisson's ratio of the material. 
Figure 2.4 shows 6 as a function of a for several values of Poisson's 
ratio. ; 

This generation of a shear Wave ctnatan ts ized in terms of 
particle velocities (Fig. 2.5). It can also be seen that a shear 
wave generated by the oblique reflection of a longitudinal wave is 
polarized in the plane of incident; that is, the particle velocity 


associated with this shear wave is normal to the direction of propaga- 


tion and directed in only one direction. The case of a compression 


le 
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Figure 2.3 dblique incidence of a plane longitudinal wave at a free 
surface 
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Figure 2.4. 
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Dependence of Angle of Reflection of Generated Shear Wave 
upon Poisson's Ratio and Angle of Incidence of Longitudinal 
Wave 
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Figure 2.5. Particle Motion Associated with an obliquely incident 
plane longitudinal wave 
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wave striking the free surface of a semi-infinite body will be 
considered, Let o be the instantaneous value of normal stress of the 
incident wave, o* the corresponding normal stress of the reflected 
longitudinal wave and t' the shear stress of the induced transverse 
wave. These stresses are all related by a coefficient of reflection R 


such that [Ref. 3] 


(2-21) o' = Ro 


(2-22) Th 


((R + 1) Cot2p Jo 


and the coefficient of reflection R is given by [Ref. 3]. 
tan 86 tan@2¢ - tan a 
(2-2.3) R= 


tan g tan-28 + tan a 


Ris a function of a and g. Since Ris a function of 8, it is actually 
a function of the Poisson's ratio of the material. Figure 2.6 shows 
R as a function of a for several values of Poisson's ratios. 

A very similar series of statements can be made for the oblique 
reflection of a shear wave. By referring back to R as defined by 


equation 2-23. 
(2-24) t' = Rt 


(2-25) of = [(R - 1)tan¢e]1 


where t, t' and o' are the instantaneous values of incident shear 
Stress, generated shear stress and generated longitudinal stress, 
respectively, 6 now representing the ancle of incidence. 

When an incident wave strikes a free surface obliquely, the 
surface will be displaced such that the angle of emergence a satisfies 


[Ref. 3, 4]. 
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Figure 2.6. Coefficient of Reflection as a Function of Angle of 
Incidence for Several Values of Poisson's Ratio 
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(2-26) a = 28 


The angle of emergence is also related to rectilinear components 


of the particle velocity by 
- -] 
(2-27) a = tan SE) 


where Vy is the sum of the particle velocity components perpendicular 
to the surface, vy is the sum of the particle velocity vectors parallel 
to the surface as seen in Fig. 2.5. Again it can be seen that the 
angle of emergence is a function of both the angle of incidence and 


Poisson's ratio. 


E. SURFACE VELOCITY REFLECTION COEFFICIENT 

The surface particle velocity for an incident compression wave is 
related to the particle velocities of the incident and reflected waves. 
For a plane elastic longitudinal wave, the particle velocity Vp at any 


time is given by the expression [Ref. 4] 


where o represents the normal stress associated with the wave. For 

a compressive wave the particle velocity is in the same direction as 
the velocity of wave propagation; for a tension wave the particle 
velocity is directed oppositely from that of propagation of tne wave. 
A corresponding relationship exists for particle velocity v, within a 


transverse wave, namely, 
(2-29) ee t/pCo 


where t is the shear stress associated with the transverse wave. 


From Fig. 2.5, the surface particle velocity is 
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(2-31) b= —— (of, + o 


where ny is the unit vector in tne direction parallel to the incident 
longitudinal wave, no is the unit vector in the direction parallel 

to the reflected longitudinal wave and na is the unit vector in the 
direction perpendicular to the reflected transverse wave. By substitut- 


ing equations 2-21 and 2-22 into the equation 2-31 


Cc 
(2-32) ¥ - = [i + AR + A, = (R + 1)Cot 28] 


or 
(2-33) Vv. = Ro 
wnere 


Cc 
(2-34) ae - [i, + A,r + fi, (R + 1)Cot 28] 


From Fig. 2.5 


ny { Cos a + j Sin a 


no -} Cos a + j Sin o 


i Sin g + j Cos B 


By substituting these equations into equation 2-34, it can be shown 


that 





a C 
(2-35) R = Hic - R) Cos a + a (1 + R) Cot 28 Sin Bj + 
] 
2 


C 
+ J[(1 + R)(Sin ao + — Cot 28 Cos if 
2 


On the other hand 


(2-36) \ = vi(i Cos 28 + j Sin 28) 


By using equations 2-33, 2-35 and 2-36, it can be shown that 


7 _ (1+ 2) Cot 28 Cos B+ Sin 8 
(2-37) R Cp (Se and 
Vp = Ro 


where R is called the surface velocity reflection coefficient. From 
this formula, it can be seen that the surface velocity reflection 
coefficient is a function of the angle of incidence, the Poisson's 


ratio and the density of the material. 
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TTI. INVESTIGATION QF SURFACE MOTION 


A. NATURE OF INVESTIGATION 

Tne theories of oblique reflection of elastic waves at a free 
Surface were outlined in tne previous chapter. Such waves cause a 
particle motion on the surface. The motion of the surface will cause 
the reflection characteristics of an acoustic pulse to change as 
various portions of the pulse encounter to surface. The theoretical 
investigations of the surface motion in solid caused by reflection 
of an intense acoustic pulse at a pressure release surface were 
facilitated by using the surface velocity reflection coefficient R 
which was derived in equation 2-5. 

The intent was to illustrate the surface motion for different 
materials, for different loadings, for different angles of incidence. 
For this purpose the IBM-360 computer was used. The surface profiles 
were obtained by using two different approaches. Tne first approach 
was the conventional one which is called the fixed surface approxima- 
tion. This assumes that surface motion can be neglected during the 
period that the pulse strikes the surface. The second approach 
is consider the motion and distortions of the surface. This is 
called the moving surface approximation. The results of these two 
approximations were tnen compared, 

The computer programs used in this investigation are listed in 
Appendix A. The fixed surface approximation computes the surface 
distortion by simply integrating the surface velocity, obtained by 
using a constant surface velocity reflection coefficient, resulting 
from a pulse arriving at the x = 0 plane. The moving surface approxima- 


tion makes two major corrections. First, the surface velocity is 
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computed at each time from the displaced position of the surface, 

the displacement resulting from the previous portions of the pulse. 

Secondly, account was taken of changes in the angle of incidence 

caused by distortions of the surface from a plane. This was done 

by constructing the normal to an equation describing the surface 

obtained by passing a quadratic curve through the point in question 

and its neighboring points. Because of tne time and position 

dependent changes in the pulse shape, numerical integration was 

used to convert the surface velocity into surface displacement. 
Throughout this investigation a saw-toothed shape stress pulse 

with 1 nm. pulse width has been used to illustrate the surface motion. 

The amplitude of the stress was changed from a value equal to the normal 

fracture stress to ten times this quantity. The normal fracture stress 

is the maximum tensile stress which the material will tolerate. A more 

exact definition would be the minimum dynamic tensional stress required 

to rupture the material. The angle of indidence was varied from 5 to 

85 degrees, using 5-degree increments. Four materials, aluminum, steel, 

plexiglass, and lead, were investigated. The characteristics of these 


materials are given in Table I. 


Beene SULTS 

Comparison of Approximations: As an illustration of the differences 
between the two approximations, the surface motion produced in aluninum 
and lead for a savtootn compressive pulse of amplitude equal to tne 
- fracture stress incident at 30° computed in the two approximations is 
illustrated in Figs. 3.1] and 3.2. The maximum displacement produced 


by the pulse is the same in the two approximations. This result is 
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Figure 3.1. Surface motions for aluminum under two different assumptions. 
Amplitude of stress is equal to normal fracture stress, a = 30°, 
time = 0.251 x 10-© sec. 
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fmrgure 3.2. Surface motions for lead under two different assumptions. 
Amplitude of stress is equal to normal fracture Stress, a = 30°, 


time = 0.741 x 107 


expected since the net displacement can be computed from the impulse 
contained within the incident pulse and the elastic moduli of the 
medium within which 1t propagates. The moving surface approximation 
predicts a more rapid displacement of the surface than the fixed 
Surface approximation. This also is quite understandable. 

The Angle ¢: The cnange of the angle of incidence due to surface 
distortion is called ¢. During the motion of surface $ at any given 
spot changes with time. To illustrate the change in » by the time at 
a given position, the curve is drawn for aluminum and given IN Figoeoeee 
It shows that the angle ¢ ataipis to increase after the pulse strikes 
to the surface, and reaches its maximum value, then decreases again. 
But the decrease is relatively slower than the increase. This is a 
result of the assumed sawtooth pulse. The angle » is also effected 
by the angle of incidence and by the amplitude of stress. The 
effect of the angle of incidence is illustrated in Fig. 3.4. 

Figure 3.4 shows the change in angle ¢ by tne angle of incidence for 
all] materials investigated. The angle » cnanges linearly with stress 
amplitude because of using the Hooke's Law and the relatively smal] 
displacements. The change in the angle $ with stress amplitude is 
given in Fig. 3.5. In this figure, aluminum is chosen as an example. 

In the moving surface assumption, the distortion of the surface 
from a plane effects the pulse shape and causes distortion on the pulse. 
Assuming plane wave propagation, there is no distortion at the front 
and tail end of the reflected pulse since it is reflectad from a flat 
surface. But in between these two ends, the different portions of the 
pulse cause different distortions on the surface which result in 


different angle of correction 9. Because of different 9, the reflected 
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Figure 3.3. 


Change in angle with time for aluminum. 
equals to normal fracture stress. 
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Figure 3.4. The change in angle ¢ with angle of incidence a for different 
materials. The amplitude of pulse is normal fracture stress. 
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Figure 3.5. 





The change in angle ¢ 


wae with stress amplitude for aluminum. 
a = 60°, t = 0.251 x 10-6 sec, 
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angle is not the same as the instant angle of incidence for these 
portions of pulse. Tnus the wave fronts for these portions of pulse 
are not in the same direction relative to the front and tail end of 
the pulse. Hence the reflected pulse under the moving surface 
assumption is not tne sane as tne incident pulse. From the investiga- 
tion, it was seen that the pulse became far apart after the reflection 
from the surface. And this causes the distortion in the stress field 
which is not investigated in here. 

The moving surface assumption also effects the amplitude of the 
reflected pulse. The amplitude of the reflected stress pulse equals 
to the product of the coefficient of reflection and the amplitude of 
the incident stress pulse. But from Fig. 2.6, it can be seen that 
the coefficient of reflection is the function of the angle of incidence. 
Since the assumption causes a change in the angle of incidence, the 
amplitude of the reflected pulse is not the same as the amplitude of 
the reflected stress pulse under the fixed surface assumption condition. 
It is either weaker or stronger than it. It depends on the angle 
of incidence and the material used. To illustrate the effect of ¢ 
on the amplitude of stress pulse, two numerical examples were done 
for aluminum. The coefficients of reflection were calculated by 
using the equation 2-23 for two different angles of incidence, 30 
and 40 degrees. The amplitude of stress was the normal fracture 
stress and the maximum angles of correction were used in both cases. 
For 30 degrees angle of indidence, by using fixed surface approximation, 
the coefficient of reflection equals -0.778. But in the same condition 


by using the moving surface approximation, it equals -0.790. For 


Sy) 


45 degrees angle of incidence, the first approximation results in the 
coefficient of reflection equals to -0.579, the second approximation 
results in the coefficient of reflection equals to -0.590. 

Surface Displacement: Surface displacement is effected by the 
Stress amplitude and by the angle of incicence. The surface displace- 
ment increases linearly with stress amplitude. However, it decreases 
by increasing the angle of incidence. The surface motions for two 
different materials, aluminum and lead, are illustrated in Figs. 3.6 


and 3.7. 


38 





ae 


ya 


Y(x107 7m) 





00 02 04 - 06 08 10 


/ 


X(x10° ‘m) 


Figure 3.6, Surface motions for aluminum. Stress amplitude equals to 
normal fracture stress and a = 30. 
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Figure 3.7. Surface motions for lead. Stress amplitude equals to the 
normal fracture Stress a = O°. 
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